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ABSTRACT 



Aims. 4U 2129+47 was discovered in the early 80's and classified as an accretion disk corona source due to its broad and partial X-ray 
eclipses. The 5.24 hr binary orbital period was inferred from the X-ray and optical light curve modulation, implying a late K or M 
spectral type companion star. The source entered a low state in 1983, during which the optical modulation disappeared and an F8 IV 
star was revealed, suggesting that 4U 2129+47 might be part of a triple system. The nature of 4U 2129+47 has since been investigated, 
but no definitive conclusion has been reached. 

Methods. Here, we present timing and spectral analyses of two XMM-Newton observations of this source, carried out in May and 
June, 2005. 

Results. We find evidence for a delay between two mid-eclipse epochs measured ~22 days apart, and we show that this delay can be 
naturally explained as being due to the orbital motion of the binary 4U 2129+47 around the center of mass of a triple system. This 
result thus provides further support in favor of the triple nature of 4U 2129+47. 

Key words, accretion: disk - binaries: eclipsing - stars: individual (4U 2129+47) -stars: neutron - X-rays: stars 



1. Introduction 



4U 2129+47 was discovered by iForman et al.1 d!978l) at a flux 
level variable between 2.4 and 4.8 x 10~ 10 erg cirT 2 s _1 (2- 
10 keV band). Observations of 4U 2129+47 in the early 80's 
showed that both its X-ray and optical light curves were mod- 
ulated over a 5.24 hr period, with a partial V-shaped min- 
imum maintaining approximately the same shape and phase 
dThorstensen et all I Il979t lUlmer et all 1 19801 IMcClintock eTal] 
fl982T hereafter MC82). A late K or M spectral type compan- 
ion of ~0.6 M was suggested, assuming it fi lled its Roche lobe, 
and t he discovery of a type I X-ray burst dGarcia & Grindlavl 
119871) led to the classification of 4U 2129+47 as a neutron star 
(NS) low mass X-ray binary ( LMXB) system dThorstensen et alj 
Il979t IMcClintock et aljfl98ll; MC82). The source distance was 
estimated to be ~ 1- 2 kpc, correspondin g to an X-ray luminosity 
of ~5xl0 34 erg s _1 dHorne et al.l lT986). The optical light curve 
could be understood in terms of the varying viewing geometry of 
the X-ray heated face of the companion, while the V-shaped min- 
imum in the X-ray light curve was interpreted as being due to the 
gradual eclipse of an extended accretion disk corona (ADC). The 
shape of the partial X-ray eclipse and the rapidity of its ingress 
and egress have been used to place upper limits on the size of 
this X-ray scattering region (~5xl0 10 cm for the 4U2129+47 
high luminosity state, MC82). The origin of ADCs is not well 
understood yet, but it is most likely related to systems in which 



the mass accretion rate is sufficiently high that a tenuous scat- 
tering corona is formed as a consequence of matter evaporation 
from the accretion disk dWhite & Holtlll982l) . 

4U 2 129 +47 was first revealed in a low state 
(Fq,3 t 6 keV^10~ 12 erg cm' 2 s~ ] ) in September 1983 dPietsch et al.l 
1986; Wenzel et al.ll 1983b . Optical observations carried out be- 
tween 1983 and 1987 showed a flat light curve without any 
evidence for orbital modulation, while the spectrum displayed 
features fully compatible w ith a late type F8 IV star dKaluznvl 
1988; IChevalier et al]|1989l) . The hypothesis of a foreground or 
a background star seemed unlikely, due to the low probability 
(<10~ 3 ) of chance superposition. This led to the 



suggestion 
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that 4 U 2 129 +47 is part of a triple system dThorstensen et al.l 
119881). The revised estimate of the source distanc e was ~6.3 kpc 
dCowlev & Schmidtkelll990l;lDeutsch etaL1ll996l) . 

Hints of a possible detection of a dynamical in teraction be- 
tween the F star and 4U 2129+47 were d iscussed by Garci aeTail 
d 19891) and lCowlev & Schmidtkel i 1990b . after the discovery of 
a ~40 km s _1 shift in the mean radial velocity measurement de- 
rived from the F star spectrum. Shifts of this amplitude are in- 
deed e xpected if the F sta r is in a month-long orbit around the 
binary dGarcia et al.lll989h . 

ROSAT and Chandra observations, carried out between 
1991 and 2000, led to a c haracterization of t he low lumi- 
nosity state of 4U 2129 +47 dGarcia et al.1 Il992t iGarcial [19941; 
iGarcia & Callananl ll999). The refined Chandra posi tion turned 
out to be coincident with the F star to within 0".l (Nowa ket al.l 
120021 hereafter N02) providing support in favor of the triple na- 
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Table 1. Mid-eclipse epoch measurements. 



Observatory 


Mid-eclipse Epoch" (JD) 


Orbital Period" (s) 


References 


L±Cjr\\J 1 dllll ZjIC/C \JU!s(Zt VUlUf y icltZSLUlsC 




1 OO J / } 


Ttinrstpnspn pt al C1 Q7Q , 1 l' r TTi7Q , l 


MrCrvnxv Mill tpJ fzcnnp anH KPNO tpJpvmnp 


2444107 785(3*1 


18857 5fn 


Mcriintock et al n981HMr8n 


Einstein 

Louisiana State Telescope, 
McGraw Hill telescope 


2444403.743(2) 


18857.48(7) 


McClintock et al. (1982) (MC82) 


Chandra 


2451879.5713(2) 


18857.631(5)* 


Nowak et al. (2002) (N02) 


XMM-Newton 


2453506.4825(3) 


18857.594(7) 


this work (tw) 


XMM-Newton 


2453528.3061(4) 




this work (tw) 



" Numbers in parentheses are the uncertainties on the last significant digit (errors at ltr level). 
b Average orbital period calculated by using the two Chandra eclipses. 



ture of 4U 2 129 +47. However, a firm conclusion could not be 
reached. 

Here we report on XMM-Newton observations of 
4U 2 129+47, and discuss the likely detection of a mid-eclipse 
epoch variation between two ~22 day distant observations. 
We show that this delay is naturally explained as being due 
to the orbital motion of the binary 4U 2129+47 with respect 
to the center of mass of a triple system. This delay is thus 
probably the first "Doppler" (or, more accurately, "Roemer") 
X-ray signature of the triple nature of 4U 2 129 +47. We outline 
our data reduction procedure in Sect. [2] and present the results 
of timing and spectral analysis in Sect. [3] Our conclusions are 
summarized in Sect. [4] 

2. Observations and data 

XMM-Newton dJansen et alj 1200 ll) observed 4U 2129+47 on 
May 15 and on June 6, 2005 for a total time span of ~80 ks 
(about four orbital periods). The total effective exposure time for 
each observation was ~ 1 3 ks for the EPIC-PN, EPIC-MOS 1 , and 
EPIC-MOS2 cameras. The remaining observing time was dis- 
carded due to ground station anomalies and high radiation from 
solar activity filling up of the EPIC-PN scientific buffer. Heavy 
contamination due to solar activity resulted in poor orbital phase 
coverage, especially during the first observation. Furthermore, 
the EPIC-PN and EPIC-MOS cameras were found to be un- 
equally affected by this contamination, thus forcing a differ- 
ent selection of good time intervals for the spectral and timing 
analyses (see Sect. |3J- The observation data files (ODFs) were 
processed to produce calibrated event lists using the standard 
XMM-Newton Science Analysis System (SAS 7.0). We used the 
epchain and emchain tasks for the EPIC-PN and the two MOS 
cameras, respectively. Source light curves and spectra were ex- 
tracted in the 0.2-10 keV band, using circles of ~14.6" radius 
centered on the source. This corresponds to -70% encircled en- 
ergy Q for both the EPIC-PN and EPIC-MOS cameras. Larger 
circles could not be used due to the proximity of the S3-/5 Digital 
Sky Survey stellar object (N02). We extracted the background 
light curves and spectra from circles of radii ~1 16" in the near- 
est source-free region to 4U 2129+47. Background and source 
circles were all chosen to lie within the same CCD. The differ- 
ence in extraction areas between source and background was ac- 
counted for by using the SAS backscale task for the spectra and 
the lcmath task from Heasoft (version 6.1.1) for the light curves. 



The average source count rate was found to be 0.041+0.001 
count s _1 in the EPIC-PN and 0.010+0.001 count s -1 in the two 
EPIC-MOS cameras (errors are ltr). Given the short effective 
exposure time, the low count rate of the EPIC-MOS 1 and EPIC- 
MOS2 cameras did not contribute significantly to the spectral 
analysis; therefore in Sect. 13.21 we discuss only the spectrum 
from the EPIC-PN camera. 



3. Results 

3. 1 . Orbital ephemerides and eclipse parameters 

In the XMM-Newton EPIC-PN light curves, two eclipses were 
clearly detected, one in each pointing. These eclipses were also 
unambiguously found in the EPIC-MOS 1 and EPIC-MOS2 data. 
In order to avoid any non-synchronicity problem s between the 
three E PIC cameras, we followed method III of Ba rnard et al.l 
d2006l) . In accordance with this method, all the source and back- 
ground time series of the same observation were extracted within 
the time interval around the eclipse that was found not to be in- 
terrupted by the presence of a previously removed solar flare 
(Sect. [2]). In all cases, the time interval selection was carried out 
by filtering each light curve with the evselect (version 3.59) key- 
words "timemin" and "timemax". This additional reduction of 
the effective exposure time (in addition to the one described in 
Sect. |2]l was especially restrictive for the first observation, for 
which a total exposure time of only ~4 ks around the eclipse 
could be used. The times of all light curves were corrected to the 
barycentre of the Solar System with the SAS barycen task (ver- 
sion 1.17.3), and summed up in each observation with the lcmath 
task, in order to maximize statistics and thus improve the accu- 
racy with which the eclipse parameters could be determined. To 
estimate the mid-eclipse times, the light curves were rebinned in 
300 s binfE 

These light curves were then fit with a rectangular eclipse 
model, in which the mean count-rate outside (F max ) and in- 
side (F m in) eclipse, and the mid-eclipse epoch (To), were treated 
as free parameters. In these fits we fixed the duration of the 
eclipse at the value 1523 s, as measured by N02. Being dic- 
tated by the secondary star, the duration of the eclipse limb is 
unlikely to have changed since the time of the Chandra ob- 
servation (see Sect. 13.21 for details), x 1 minimization was per- 
formed with an IDL routine written by the authors. The model 
rectangular eclipse was integrated over each time bin before 



1 As described in chapter 3.2.1.1 (Issue 2.5) of the XMM-Newton 
users' handbook. 



2 A check was carried out a posteriori to verify that our results are 
virtually independent of the light binning time. 
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Fig. 1. Fit of the mid-eclipse epoch during the first observation. 
The 0.2-10 keV light curve (bin time 300 s) is shown together 
with the best fit model (dashed line). 



Fig. 2. Fit of the mid-eclipse epoch during the second obser- 
vation. The 0.2-10 keV light curve (bin time 300 s) is shown 
together with the best fit model (dashed line). 



the x 1 was computed, in order to take data binning into ac- 
count With this method we obtained an accurate determi- 
nation of the eclipse mid-epoch, even though the ingress and 
egress eclipse times could not be determined with the same 
accuracy. The x 2 hyper-surface was directly sampled in a fine 
grid of values in order to distinguish local minima. The vari- 
ance between model and data was then calculated in each point 
and for each set of parameters, in order to investigate the lo- 
cal x 2 minima in the 4D parameter space. The best fits to the 
eclipse epochs were found to be T (a)=2453506.4825+0.0003 
JD and T (b)=2453528.3061±0.0004 JD with ;f 2 /dof(a)=9/l 1, 
^ 2 /dof(b) =44/42 (errors are at 67% confidence level unless oth- 
erwise specified; our epochs are given in UT0). The values of 
the reduced x 2 in the above fits are close to 1 and therefore the 
addition of any other free parameter in the fit would not be jus- 
tified from a statistical point of view. We also checked that by 
allowing the eclipse duration to vary within the N02's lcr con- 
fidence level, i.e., 1473-1553 s, the other parameters of the best 
fit remain unchanged to within the errors (mid-eclipse epochs 
differed by less than 10 s and the lcr errors remained below 26 
s). We note however that if the eclipse duration is included as a 
free parameter without any constrains, a less significant eclipse 
delay of <5r m =151+55 s would be obtained. 

In the following we adopt values of the two mid-eclipse 
epochs as derived from the first fit. We show in Figs. Q] and 
[2] the two eclipses, together with the best fit models discussed 
above. We also carried out the fits by u s ing a modified version 
of the eclipse model by iMangano et al.l d2004l) in order to esti- 
mate the eclipse ingress and egress time. We allowed these times 
to take independent values during the fit (N02). Only upper lim- 



3 Standard fitting routines that compute the fitting function punctu- 
ally in the center of the bin are not adequate when a function with a 
large first derivative or features with a scale smaller than the bin time is 
considered. In our case the derivative diverges at the ingress and egress 
times. 

4 Note that HJD(UT)=HID(TT)-64.68 s at our epochs. We did not 
consider the correction for the difference between heliocentric lulian 
date in the geocentric (terrestrial) dynamical time system, HID(TT), 
and barycentric dynamical time system, BID(TB). The latter is the one 
used by the SAS barycen task, but the difference between BID(TB) and 
HID(TT) is less than ~3 s at any given time, which is much smaller than 
the accuracy of our measurement here. 



its of <3 10 s could be derived, which are significantly larger 
than those obtained from the Chandra observation (in 130-260 
s range, N02). 

In order to determine a refined orbital solution, we consid- 
ered the above mid-eclipse epochs together with the epochs, T n , 
derived from earlier observations (see Table [T). As discussed in 
Sect. [T| the observed phase alignment between the X-ray and 
optical light curve minima, allows the comparison of optical 
and X-ray measurements of the system's ephemerides. The long 
time span covered by the eclipse measurements (1979-2005 or 
«max~44800) can be used to improve the accuracy of the orbital 
solution and, possibly, measure the orbital period derivative. To 
this aim we used a standard O-C techniquqfl We considered the 
ephemeris from MC82 as reference (r ref =2444403.743+0.002 
JD, P ret -= 18857.48+0.07 s), and plotted in Fig. E the delays 
AT n =T n -T u ^. Here r nprcd =r re f+nP re f, with n the closest integer 
to (T n -T re f)/P le f (our two observations correspond to n=41706, 
41806). In the same figure we have also plotted the best 
quadratic fit to the O-C residuals, corresponding to an orbital 
period evolution with constant time derivative. A linear fit (i.e., 
a constant orbital period) to the same data gave an unacceptable 
fit (^ 2 /dof=91/4). Table [2] gives our corrected reference time, 
r re f, orbital period, P re f , and the derived orbital period evolution 
P orb p- r J j =(5.8+0.7)xl0 6 yr. This value is a factor of ~4 larger 
than that in N02, but we note that N02's estimate was deduced 
by adopting the MC82 value of P or b and accounting for the entire 
measured delay as being due to an orbital period derivative. The 
observed delay of ~6500 s between the ephemeris of MC82 and 
the one found in the present work (see Fig. [3]), implies an orbital 
phase shift of ~0.35. Small deviations in the eclipse centroid 
between the active and quiescent state of 4U 2129+47 may 
be introduced by changes in the shape of the ADC around the 
compact object. However a ~6500 s delay is far too large to be 
explained as a result of such geometrical variations. The poor^ 2 
in Table|2]is due to the large shift of the first XMM-Newton point 
with respect to the second one (£7' m =192±43 s), which could 
not be accounted for by any quadratic fit to the eclipse phase 
evolution. While this delay is much smaller than that discussed 



5 Observed minus calculated residuals, which are the delay in 
eclipse time over that expected for a c onstant period system (see, e.g., 
Parmar et al. 1991 ; Papitto et al. 2005) and references therein. 
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Fig. 3. Delays of the mid-eclipse epochs with respect to a con- 
stant Porb— -fW model. The solid line in the upper panel repre- 
sents the quadratic best fit to the epochs (see Sect. 13. It . The 
lower panel shows the residuals from this fit. 



above, we argue that it is also very unlikely to result from geo- 
metrical variations within the 4U 2129+47 binary system. This 
is because the eclipse profile is consistent with the central source 
being eclipsed by the companion star. Moreover, the source 
X-ray flux and spectrum remained virtually the same across the 
two XMM-Newton observations and the Chandra observation 
discussed by N02 (see Sect. I3.21 >. In Sect.|4]we discuss the possi- 
bility that this delay is due to light propagation in a triple system. 

We also extracted the EPIC-PN light curves of the ob- 
servations, since these have a better orbital phase coverage 
than those obtained by summing all three instruments. This 
however resulted in a lower count rate and S/N. These light 
curves and those obtained by using data summed over the three 
EPIC cameras were folded at the best orbital solution using 10 
phase bins. We fitted these two light curves with the function 
F(0)=A+Bsin[27r(0-0o)], and looked for a sinusoidal modulation 
similar to that observed by N02. No significant modulation was 
observed (fitting with a constant value gave a^- 2 /dof of 12/13 and 
10/17 respectively). We derived a 90% confidence upper limit on 
the amplitude modulation of ~17%, i.e., a factor of ~2 smaller 
than the value reported by N02. This result was also checked by 
using light curves extracted only in the 0.5-2.0 keV band, where 
the amplitude modulation might be higher (N02). No significant 
differences were found. We discuss these results in Sect. [4] to- 
gether with the results from the spectral analysis (Sect. [3~!2l . 

3.2. Spectral analysis 

Spectral analysis was carried out by using XSPEC version 
11.3.2t (Arnaud 1996); the data were rebinned in order to have 



Table 2. Orbital solution obtained with the best quadratic fit to 
the O-C delays AT n (see Fig. 0. 



(JD) 


2444403.7443+0.0013 


forb (S) 


18857.594+0.007 


Porb (S S-') 


(1.03±0.13)xl0-'° 


PorbP^ (yr) 


(5.8±0.7)xl0 6 


X 2 /dof 


25.6/3 



o o 




it 



0.5 



Energy (keV) 



Fig. 4. Measured 0.2-10 keV spectrum of 4U 2129+47 during 
the second pointing. The best fit model and the contribution of 
the fit residuals to the x 2 are also shown. 



at least 20 photons in each energy bin. Owing to poor statistics, 
phase-resolved spectroscopy could not be carried out. Instead, 
the spectra of the two observations were accumulated during the 
same time intervals selected for the extraction of the EPIC-PN 
light curves, except for the eclipses (see Sect. 13. U . 

The spectrum of the first observation was modeled with only 
an absorbed blackbody (poor statistics did not allow for more 
complex models). 

The second observation was first modeled with an absorbed 
power law plus a blackbody component, but the F-statistics prob- 
ability with respect to a simple absorbed blackbody model was 
found to be ~0.2. The best fit was then obtained, also in this case, 
by adopting an absorbed blackbody model. 
A power law component with fixed photon index T= 1 . 1 (N02) 
was added to the fit of the second observation in order to esti- 
mate an upper limit. For such a power law component, the 90% 
confidence upper limit was about 10% of the 0.2-10 keV un- 
absorbed flux (in agree ment with N 02, see their model A). C- 
statistics model fitting dCashl [19791) was also performed on the 
unbinned spectra: the results were fully compatible with those 
obtained by using x 1 minimization. Figure [4] shows the spec- 
trum and model of the second pointing as an example, while the 
best fit parameters are reported in Table [3] No significant dif- 
ference was found between the two XMM-Newton observations, 
and all parameters were compatible, to within the errors, with 
those determined by N02 for the quiescent state of 4U 2 129+47. 
We a l so tested our resu lts with the xspec nsa model dArnaudl 
119961; IZavlin et al.lll996l) . Fits were carried out, first by using 
a fixed distance of 6.3 kpc, and then by fixing a neutron star ra- 
dius of 5 and 10 km (N02). Results of these fits were found to 
agree with those of N02 (our errors on all parameters are a factor 
of -1.5 larger). 

4. Discussion 

We reported on XMM-Newton observations of 4U 2129+47 in its 
quiesc ent state, which has lasted, apparently uninterrupted, since 
1983 (IWenzel et al.lll983l) . The discovery of a late F-type star 
coincident with the positio n of 4U 2129+47 dThorstensen et al.l 
ll988HChevalier et al.ll 19891) led to the hypothesis that this binary 
system might be part of a hierarchical triple. Our detection of a 
delay 6T m - 192+43 s across two eclipses separated by ~22 days, 
can be naturally explained as being due to the orbital motion of 
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1 Obs. 


(May- 15-2005) 


2 Obs. 


(June-6-2005) 




* 2 -Stat 


C-Stat 


* 2 -Stat 


C-Stat 


JVh (10 22 cm" 2 ) 
kT bb (keV) 

K b ( km ) 
X 2 /dof 

C-Stat 

^-lo^dO-'^rgcm- 2 s- 1 ) 
^Lio ke v(10- 13 ergcm- 2 s" 1 ) 
^o.2-io k eV do 32 ergs-') 


0.14 + °'i 

-0.09 

25 +0M 

14.30/14 

9.3 
1.6 

7.7 


0.24 + »°3 

" -0.09 

22 +0 03 

"•■"-002 
1 Q+0.74 
1 - 7 -0.53 

496.8 
9.3 
2.3 
10.9 


0.21 + °L 

-0.06 

20 +0 02 

-0.02 
7 4.7+1.8 
^■•^-0.7 

20.6/23 

8.9 
2.2 
10.5 


-0.(8 

2.87^ 

592.24 
9.0 
3.3 
15.8 



a Absorbed Flux 
b Unabsorbed Flux 

c From the unabsorbed flux and assuming a distance of 6.3 kpc. 



the binary with respect to the center of mass of a triple, and thus 
lends support in favor of the triple system hypothesis. 

Using a third star of mass Mi, an inner binary with M 2 ~2 
M , and a non-eccentric orbit, the expected delay between two 
mid-eclipse epochs separated by a time interval t can be ex- 
pressed as 

5T = a-i sin ij c [sin (</>o + 5<p) - sin <po] . ( 1 ) 

Here a 2 =(G/4^ 2 ) 1/3 P 2/3 Mi/(M l +M 2 ) 2/3 is the radius of the in- 
ner binary orbit around the triple system center of mass (with a 
period P tr ), </>o is the phase at 7o(a) of the binary system along 
such an orbit, c is the speed of light, i is the inclination angle of 
the triple system orbit, 5<p=2nT and, in our case, r=7o(b)- 
T (a)^22 d. 

Figure [5] shows a plot of the range of allowed values of 
the triple orbital period as a function of the third star mas^l, 
which give delays compatible with the measured value 6T m (to 
within the uncertainties at a given confidence level). The dashed, 
dot-dashed and dot-dot-dot-dashed lines represent the ranges for 
l<x=43 s, 2<x=85 s and 3<x=165 s confidence intervals in 6T m , re- 
spectively. The solid lines give the constraints on P tr imposed by 
the me asured ~40 km s -1 shift in the mean radial velocity of the 
F star dGarcia et alJI 19891) . under the assumption that this repre- 
sents the maximum observable radial velocity shift of its orbit. 
The two panels of Fig. [5] show the cases of i=90°, and i=60° re- 
spectively. From the upper panel of this figure (i=90°) it can be 
seen that, at Icr confidence level, we can set a lower limit on the 
F star mass of ~1.2 M . Considering a 2cr uncertainty on 6T m 
removes the upper limit on Mi and only the orbital period can 
be constrained (39 d<f t r<200 d, for M^l M ). We note that 
a decrease of the inclination angle results in a smaller range of 
allowed orbital periods, while the effect of a non-zero eccentric- 
ity (hypothesis not considered in our calculation) would have the 
opposite effect. 

The possibility that 4U 2129+47 is in a triple system might 
also have noticeable consequences for our measured orbital 
period derivative, P Qrh p-^=(5.S±0J)xl0 6 yr (see Sect. UTi 
At first sight, this value se ems to imply that, as in other known 
LMXBs (e.g., XI 822-37 1 , iParmar et alj|2000t iHeinz & Nowakl 
l200lt IWhiteetalJ[l995l) the orbital period of 4U 2129+47 is 
increasing. This would be contrary to evolutionary expectations. 
In fact, angular momentum losses (such as gravitational wave 

6 For an F-type star it is expected Mi~l-1.6 M (see, e.g., 
iBohm-Vitensell 19921) . 



emission and magnetic breaking) in a binary system with an 
orbital period of ~5 hr and a non- degenerate com panion would 
imply a decreasing orbital period dVerbuntl ll993). However, if a 
modulation in the eclipse arrival time of amplitude r a =«2sin i/c 
(as discussed above) contributes to the measured epochs of the 
eclipses, caution should be used with any estimate of P or hP^ h 
(such as that in Sect. 13. U that does not take this modulation 
into account. Given our poor knowledge of the parameters of 
the triple system, a detailed correction for the orbital motion of 
the inner binary around the center of mass of the triple cannot 
be carried out yet. A conservative estimate of P or bP^ h can be 
derived by increasing the uncertainties on the observed AT n 
up to a value T a , which represents the unknown amplitude 
of the modulation in the eclipse arrival times. The range of 
values that T a may attain depends mainly on the range of 
orbital periods for which a solution of EqJTJ compatible with 
the values of r and 6T m we measured, exists. Unfortunately, 
the range on P tr that can be inferred from Fig. [3] is fairly 
loose. Based on our Icr range of 6T m , we derive a minimum 
upper limit of P" pp ^75 d for the third star orbital period (for 
Mi=d.2 M and i=90°), which corresponds to r a ^100 s. 
We checked that increasing the uncertainty on the measured 
values of the eclipse epochs by this value would remove 
entirely the need to introduce an orbital period derivative. By 
analogy we speculate that the presence of a third body in a 
wide orbit around the inner binary might explain the positive 
orbital period derivatives observed in some LMXBs, which is at 
odds with the expectations of the standard evolutionary scenario. 

Due to the low count-rate, and short effective exposure time, 
our XMM-Newton observations did not detect any spectral com- 
ponent above 2 keV. An absorbed (Af H ~0.2xl0 22 ctrT 2 ) black- 
body component with kTbb~0.2 keV and /?bt>~2 km provided an 
adequate modeling of the spectrum. 

For quiescent NSs in LMXBs, like 4U 2129+47, this soft X- 
ray emission can be produced in several alternative ways. One 
possibility is that this emission is powered by thermal energy re- 
leased as the NS cools in between accre tion phases dBrown et al.l 
119981: IColpi et alJl200lUWiinandsll2002l) . Alternative models in- 
voke a NS in the propeller regime, shock emission due to the 
interaction between the pulsar wind and matter in the vicinity 
of the compan ion star, and residual accretion onto the compact 
star (see, e.g.. [Stella et alj|1994l ICampana et alj|1998l) . While 
the properties of the propeller regime are still rather uncertain, a 
power law emission of photon index 1-2 should be expected at 
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Fig. 5. Allowed regions of the triple orbital period, P tl , as a func- 
tion of the third star mass, Mi, in the cases in which the mea- 
sured delay 5T m =192 s is varied within l<x (dashed line), 2<x 
(dot-dashed line) and 3<x (dot-dot-dot-dashed line) confidence 
level. For the 3cr range only the lower limit is drawn. Solid 
lines represent the limits imposed by the measured radial veloc- 
ity shifts of the F-star (see text for more details). The two panels 
show the cases i=90° and 60°. 



least in the last two interpretations (Campa na et al.| [l998; N02). 
Discriminating between these models is not possible based on 
our results, given the lack of information at energies above 2 
keV. 

A soft spectral component with similar properties, plus a 
power law component with photon index -1.1, was observed 
in the quiescent state of 4U 2 129 +47 in December 2000 (N02). 
However, also this observation was hampered by a low number 
of counts above 2 keV, thus preventing an accurate characteri- 
zation of the power law emission. Some indication was found 
that the power law component was consistent with being of con- 
stant amplitude and slope, while the blackbody component was 
sinusoidally modulated over the orbital period, in a manner con- 
sistent with neutral column density variation (a factor of ~2). 
This modulation was ascribed to the presence of a vertically ex- 
tended disk atmosphere, thicker at the outer rim, close to the 
region where the accretion stream from the secondary star im- 
pacts. 

As discussed in Sect. 13.21 our folded light curve showed 
no clear indication of a sinusoidal modulation. The upper limit 
we derived on the amplitude of this modulation is significantly 
lower than that discussed by N02 (a factor of ~2). This suggests 
that there has been some change in the geometry of the outer 
disk region (in particular the region where the stream from the 
secondary impacts) across the Chandra and XMM-Newton ob- 
servations. 

A series of monthly spaced XMM-Newton or Chandra obser- 
vations would afford a much more accurate characterization of 
the quiescent emission of 4U 2129+47 and measure with good 
accuracy the modulation in the X-ray eclipse times due to the 
third body. This will yield much needed information on the triple 



system parameters and single out unambiguously any orbital pe- 
riod evolution. 
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